Theory for high spin systems with orbital degeneracy 



Shun-Qing Shcni^^^ x. C. Xie^'^, and F. C. Zhang^ 
^Department of Physics, The University of Hong Kong, Pokfulam, Hong Kong 
^Department of Physics, Oklahoma State University, Stillwater, OK74-078 
^Department of Physics, University of Cincinnati, Cincinnati, Ohio 45221 
Institute of Physics, Chinese Academy of Sciences, Beijing 100080, China 
'^International Center for Quantum Structure, Chinese Academy of Sciences, Beijing, China 

(Nov. 9, 2001) 

High-spin systems with orbital degeneracy are studied in the large spin limit. In the absence of 
Hund's coupling, the classical spin model is mapped onto disconnected orbital systems with spins 
up and down, respectively. The ground state of the isotropic model is an orbital valence bond 
state where each bond is an orbital singlet with parallel spins, and neighbouring bonds interact 
antiferromagnetically. The possible relevance to the transition metal oxides are discussed. 

PACS numbers: 75.10.-b 



In many transition metal oxides, there is an orbital 
degeneracy in the electron occupation energy, resulting 
in rich and novel magnetic phenomena. |Q| The orbital 
ordering and orbital density wave have been observed 
experimentally in a family of manganites. |^ In these 
systems, the spin coupling depends on the electron's or- 
bital occupations. The simplest model to describe spin- 
i systems with two-fold orbital degeneracy is the SU(4) 
model where the Hund's rule coupling and the orbital 
anisotropy are neglected. Some transition metal 

oxides, such as manganites and lanthanum vanadium ox- 
ides, have higher spins, which has attracted much the- 
oretical interests. 1^-^ Mathematical models for high 
spin systems are generally complicated. The Hund's 
rule coupling favors ferromagnetic (FM) spins, and the 
anisotropy of the electron hopping integrals breaks or- 
bital SU(2) symmetry. While the phase diagrams of these 
complex spin systems are rich, there is lack of a tractable 
method to systematically study these systems. In this pa- 
per, we study a class of spin systems with two-fold orbital 
degeneracy in the limits of large spin. Our model can be 
derived from a lattice of ions with high spins fcrromag- 
netically coupled to the strongly correlated electrons with 
double orbital degeneracy. In the classical limit of spin, 
such a system can be mapped onto two decoupled black 
and white sub-systems, representing classical spin up and 
down, respectively. This provides an efficient method to 
study the ground states. [|l0| We apply this method to 
study the case where the electron hopping integrals are 
isotropic and use the linearized spin wave theory to ex- 
amine the effects of the quantum spin fluctuation. The 
ground state in ID chain is a spin disordered valence 
bond state, in which each bond is an orbital valence bond 
(orbital singlet, parallel spins) and neighboring bonds in- 
teract antiferromagnetically. The ground states in a 2D 
square lattice and in a 3D cubic lattice are spin antiferro- 
magnetically ordered orbital valence bond (OVB) state. 
We discuss possible relevance of the theory to some tran- 
sition metal oxides. 



We start with an Hamiltonian describing correlated 
electrons coupled to the spins of local ions. 



In the above equations, a — 1,2 are the orbital indices 
and (T is the electron spin. The sum over (ij) runs all 
the nearest neighbor (n.n.) pairs. The electronic part of 
the Hamiltonian is a generalized Hubbard Hamiltonian 
with two-fold orbital degeneracy, where the Hund's rule 
coupling and anisotropy of the electron hopping integrals 
are neglected, s'f^ and sl*""-* are the spin operators of 
electron and ion at site i, respectively. We consider Jh > 
0, consistent with the Hund's rule. In the strong coupling 
limit, U ^ t, sJh, with s the quantum number of the ion 
spin, the projection perturbation theory may be applied 
to study low energy physics of the system. At the filling 
of one electron per site, and up to the order of t'^/U, this 
leads to an effective Hamiltonian, jl^ 
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where is the electron orbital operator, J = At'^/U. 
Furthermore, in the case sJh > f^/U, the electron and 
ion spins strongly bind ferromagnetically to form a state 
with total spin st = s -I- 1/2. The effective Hamiltonian 
within this Hilbert space of total spin st can be obtained 
by applying the projection operator, p2[ 
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The new coupled spin-orbital model is found to be , 
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with Si — S^f* + s^'""-' the total spin operator with max- 
imal eigenvalue Sf at site i. 

The model has SU(2) xSU(2) symmetry, representing 
rotational invariance in both spin and orbital spaces. A 
special case is at Sf — 1/2, corresponding to s|*°"^ ~ 0. 
In that case, the Hamiltonian possesses a higher SU(4) 
symmetry and the model has been studied extensively. In 
this paper, we focus on another limit where Sf 3> 1. We 
start with the classical spin to replace Si by a classical 
vector of length s* described by two angles di and ^i. 
The Hamiltonian then reads. 
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where is the angle between the two spin vectors. To 
investigate the ground state, we use the variational prin- 
ciple to find the equations for 9i and (pi. dHc/60i = and 



dHc/5(j)i ~ 0. We see that 



and 9i 



9o or /and 



9o+n at all the sites are solutions of the equations. Below 
we will consider these solutions. In ID it can be shown 
that these solutions give the lowest energy [|l^. In 2D 
or 3D we speculate these solutions contain the lowest en- 
ergy states, although a rigorous proof is absent. We set 
9f) ^ below, and label all the lattice sites with 9i = 
by blacks, and all the sites with 9i = tt hy whites. The 
bond Hamiltonian is then reduced to 



J(T, -T, +1/4), if 0, 
0, 



if 9, = 9, ± TT. 
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Therefore, the coupled spin-orbital system is decomposed 
into disconnected black and white sites or blocks (collec- 
tion of the connected same colored sites). In the same 
colored block, all the spins are parallel but the orbitals 
interact antiferromagnetically; and the interaction van- 
ishes between different blocks. The total energy of the 
system is then the simple sum of the these blocks. This 
greatly simplifies the calculations. 

Let us first consider a two-site problem. The ground 
state of Eq.(^) is an orbital singlet with total spin 2s t- 
We shall call this 2-site state as an orbital valence bond 
(OVB), whose energy is —0.5 J. There are two compet- 
ing terms in Eq. (|^) . One is the spin-orbital coupled term, 
cos^(Oij/2)Ti • Tj, which favors FM spins and antifer- 
romagnetic (AFM) orbitals. The other is the spin in- 
teraction term cos^(Oij/2)/4, which favors AFM spins. 
Because of this competition, the ground state of more 
than 2-sites is generally not a uniform FM spin state as 
we will see explicitly below. In the ID chain, we divide 
the chain into blacks and whites and calculate the lowest 
energies of these spin configurations. In particular, we 



consider the configurations with alternating black and 
white segments of n-sites (see Fig. la). The total energy 
of the chain is the sum of these independent segments. 
The ground state energies of these segments are calcu- 
lated by using exact numerical diagonalization method 
from n = 2 to 16. The results are plotted in Fig. 2. The 
energies of the even n segments are lower than those of 
the odd n ones. The ground state is an OVB solid with 
alternating spins up and down with an energy — 0.25J 
per bond, much lower than that of the uniform FM spin 
state corresponding to n — > oo. The energy of the lat- 
ter state can be deduced from the result of the Bethe 
ansatz solution, and it is (1 — 2 In 2) J/2 = —0.193 J per 
bond. Our numerical results approach to the exact re- 
sult rapidly as n increases. Including the spin quantum 
fluctuation, we expect that the ID AFM spin long range 
order be destroyed based on the dimensionality consid- 
eration and the calculation of the linearized spin wave 
theory. The system is then described by a spin-2st chain 
of the OVBs with AFM Heisenberg couphng J/{32s^). 
Since 2st is an integer, the ground state is a Haldane's 
gap state. The quantum ground state has two-fold de- 
generacy corresponding to the translational invariance by 
one lattice constant in the OVBs. 

In a 2D square lattice, we consider various spin 
configurations including alternating colored 2-site bond 
state (Fig. lb), alternating 4-site plaquette state (Fig.lc), 
stripes and the uniform FM states. The ground states 
are found to be highly degenerate. The lowest en- 
ergy state in Fig. lb is an OVB state, and the lowest 
energy state in Fig. Ic is a plaquette orbital singlet state 
whose orbital is given by ^ t(12)t(34)-t(14)t(23), with 
t(12) representing an orbital singlet of sites 1 and 2. 
The OVB, the plaquette orbital, and their mixed states 
are degenerate, and have energy — 0.125J per bond. 
The uniform FM spin state has a much high energy 
< Ti ■ Tj + 1/4: >2D= -0.085 J per bond as estimated 
from the known result of 2D spin-1/2 model The 
OVB state is also found to have the lowest energy in 3D 
cubic lattice. The large degeneracy of the ground states 
found in the classical spin limit in higher dimension is 
removed when the spin quantum fiuctuation is included. 
We have used the spin wave theory to calculate the en- 
ergy correction to the classical spin states in 2D square 
lattice, and found it to be — 0.01157J/st for the OVB 
state, and —0.01012 J/sf for the plaquette orbital state. 
Therefore the quantum spin fiuctuation favors the OVB 
state. In 2D and 3D, the AFM spin long range order is 
expected to survive from quantum fluctuation. We thus 
conclude that the ground state of Eq.(|^) in 2D and 3D 
is the OVB state with AFM ordered spins between the 
n.n. bonds. 

We now turn to discuss the possible relevance of our 
theory to some transition metal oxides. We first exam- 
ine the ground state of cubic vanadate LaVOa. The ob- 
served magnetic order of the vanadate is of C-type AFM 



2 



phase (FM along the c-axis and AFM in the a-b plane), 
arising from a practically undistorted structure above 
the Neel temperatures. fl^Jl^ To explain the unusual 
magnetic ordering, KhaliuUin et al. ^ recently proposed 
a spin-orbital Hamiltonian for the oxides. In the limit 
rj — Jh/U — + (the notations are the same as in their 
paper) and along a given cubic axis, their model is re- 
duced to Eq.(^ in the present paper with st = 1- The 
corresponding two equivalent orbitals for the coupling 
along c-axis are, dzx and dyz, for example. KhaliuUin et 
al. compared the ground state energies of the C-type and 
G-type (AFM in all three directions) AFM phases. They 
concluded that the C-type AFM is more stable at 77 = 0, 
and is further stabilized at 77 > 0. In that work, the com- 
parison between the C-type AFM and the OVB state 
was not included. As we discussed earlier, the ground 
state is an OVB state at r; = in the large spin limit. 
We believe the experimentally observed C-type AFM in 
that compound is stabilized by the Hund's coupling. In- 
deed we have compared the energies of C-type AFM and 
OVB state in their model (Eq.(l) of Ref. [|), and found 
that the C-type AFM has lower energy for 77 > 7]c with 
r]c ~ 0.06. It will be interesting to examine the higher 
order spin fluctuation or to use numerical techniques to 
verify if the large spin limit applies to spin-1 system. 

In the second example, we consider the spin-orbital 
model for LaMnOs assuming the lattice is not distorted. 
This type of model has been studied by many authors 
and the ground state is A-type AFM (AFM along c-axis 
and FM within a-b planes). Here we wish to point 
out that the A-type ordering may be obtained by treating 
the Hund's rule coupling as a small perturbations. The 
effective Hamiltonian describing spin S — 2 and the two 
degenerate e^ orbitals of the Mn-ion is given in the limit 

A - 0, if 

where 7 is along j — i and the orbital operators = 
cos(2m^7r/3)T^ - sin(2TO^7r/3)Tf (m^ = 1,2,3) with 
eigenvalues ±i. The ground state of this Hamilto- 
nian in the large spin limit is highly degenerate. As 
the expectation value of (r^ ~^ h) {'''] + I) always not 
less than zero, the low bound for the ground state en- 
ergy is zero. Thus the G-, C-, and A-type AFM are 
all degenerate ground states with zero energy. When 
a small Hund's rule coupling is introduced, the high 
degeneracy in the ground state is removed. To sim- 
plify our problem, we keep the perturbation term due 
to the Hund's rule coupling up to the order JJ^/U, 
AH = -{J/2sl) J2 (Sj • + si) /i, where h is of order 
sJh/U < 1. Including AiJ, the A-type AFM spin state 
is most favorable in energy. In that state, the orbitals are 
C-type antiferro-orbitally ordered: dj,2_22/dy2_^2 alter- 
native in the a-b plane, and ferro-orbital along the c-axis. 



This state is favored because of the energy gain of AiJ 
from the bonds within the a-b planes. Our calculation 
shows that the superexchange interaction including the 
Hund's coupling leads to a spin A-type and orbital C- 
type structure. This state might be realized in KCuFs. 

Note that the orbital ordering pattern is different 
from that experimentally observed in LaMnOs with a 
d32;2_r2/d3y2_r2 alternating ordering. The observed or- 
bital ordering is due to strong Jahn- Teller distortion , 
but not due to the superexchange 

Finally, we examine an extended model of Eq. (^) in a 
cubic lattice to include the anisotropic coupling in differ- 
ent axis and the Hund's rule coupling, 

i/ = ^J„(S, •S,+a + s?)( T,-T, + |) (8) 

where is the coupling strength along the direction a = 
a,b,c. At the symmetric point B = 1 and Ja = Jb = Jc, 
the model is reduced to Eq.(|^) . In the presence of Jh, 
B < 1 (see in AH in the previous example). We com- 
pare the ground state energies of various states including 
OVB state, the uniform FM spin state, and the A-, C-, 
and G-type AFM. In Fig. 3, we plot the phase diagram 
for the lowest energy states in the parameter space B and 
1] — Jab/ Jc- The OVB state is found to be stable ground 
state in a finite parameter space around the symmetry 
point. It will be interesting to find an experimental real- 
ization of the OVB state. Since the OVB state break the 
translational symmetry, we expect the electron-lattice in- 
teraction may accompany a crystal structure phase tran- 
sition to result in shorter and longer bonds. At a larger 
77 and smaller S, C-type AFM spin state is realized, in 
which the orbitals are disordered liquid. At a smaller 
77 and smaller S, A-type AFM becomes stable, in which 
the orbitals has a long range AFM order in the x-y plane, 
and FM along the z-axis. The phase diagram is thus very 
rich, reminiscent of some of the features in the transition 
metal oxides. 

In conclusion, we have examined the ground state of 
a coupled high-spin orbital model from the large spin 
limit point of view. A novel orbital valence bond state is 
found to be stable in a certain parameter region, and may 
be realized in future experiments. The anisotropic model 
contains both C-type and A-type AFM spin states among 
others. Our method should be useful in study of more 
realistic models when Hund' rule coupling is weak. The 
large spin approach emphasizes orbital quantum fluctu- 
ation over spins, and special care is needed to apply the 
method to spin-1/2 case. In the SU(4) limit, there is 
a permutation symmetry between spin and orbital, and 
we do not expect the method will work. The quantum 
fluctuation may drive the two degenerated spin and or- 
bital valence bond states to a quantum state with more 
intricate correlations. In fact, from the Bethe ansatz so- 
lution, the ID chain has a gapless SU(4) liquid ground 
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state, which is not dimerized. |2^] 
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FIG. 1. Illustration of several possible lowest energy states 
for the model. The solid line represents OVB with spin up, 
and dashed line represents OVB with spin down. (a). A 
one- dimensional OVB solid, (b). A two-dimensional OVB 
solid, (c) A two-dimensional plaquette RVB solid. 

FIG. 2. Energy per bond (Js?/2) versus the site number 
per cluster on a one-dimensional chain. The solid line corre- 
sponds to the exact energy of one-dimensional chain by means 
of Bethe ansatz. 
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FIG. 3. The phase diagram for high spin systems with or- 
bital degeneracy in Eq. (5). OVB(ab) represents the OVBs 
aligning along the a-or b-axis, and OVB(c) represents the 
OVBs aligning along the c-axis. 
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